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Abstract 

^ ■ We investigate the dynamics of a four-photon Jaynes-Cummings model for large 

photon number. It is shown that at certain times the cavity field is in a pure state 
. which is a superposition of two Kerr states, analogous to the Schrodinger cat state 

(superposition of two coherent states) which occurs in the one and two photon cases. 



1 Introduction 

The Jaynes-Cummings model (JCM) is a standard and important model which describes 
the interaction of an atom and a single mode radiation field in a cavity (l], 0]. The 



importance of this model is due not only to its exact solvability, but also arises from 

P-i. 

some of its purely quantum effects, such as the periodic collapse and revival of the atomic 

"S ' 

number population, and the behaviour of the cavity field which becomes a Schrodinger 

d ! 

cat state - a phase-correlated superposition of two coherent states - at certain times. The 
JCM model has been generalized in many different ways. One such is the multiphoton 
' generalization described by the following Hamiltonian (h = 1) 



H = ua)a + -LJ a 3 + g(a) k O- + a k a + ), (1.1) 

where A; is a positive integer, 03 = |e)(e| — cr+ = |e)(<?| and er_ = \g)(e\ are atomic 

operators, g is the coupling constant, luq and u> are the atomic transition frequency and 
cavity resonant mode frequency respectively. It is interesting that for both the one jj| 
and two ||] photon cases, the cavity field is in a Schrodinger cat state at intervals 
corresponding to one half of the revival time. This arises from the fact that the Rabi 
frequency in both cases is proportional to n (up to a constant) in the approximation of 
large photon number || |], ||. 
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In this paper we shall address the dynamics and quantum characteristics of the 4- 
photon JCM in the large photon number approximation. Under this approximation the 
generalized Rabi frequency (Eq.(2.3)) depends on n nonlinearly (up to n 2 ) and this non- 
linearity causes Kerr nonlinearity in the cavity. We shall study the photon number distri- 
bution (Sec. 3) and the entropy of the cavity field (Sec. 4) and the Q-function (Sec. 5). 
It is found that the cavity field is in a Kerr state or a macroscopic superposition of two 
Kerr states at certain times. We also study the atomic number inversion and find that it 
exhibits collapse and revival for short time intervals. 

2 Large photon number approximation 

We assume that at the initial time t = 0, the atom and field are decoupled and the atom 
is initially prepared in the excited state |e), while the field is in the coherent state \a) 



CO 

|M)=^C n |n), C„ = e ~' tt ' 2 /2« (2.1) 



n=0 

where a is a complex number. For simplicity, we only consider the on-resonance interaction 
case Wq = kuj as in B . Then the combined atom- field wave function at time t is obtained 

as 

oo 

l*( T )> = J2 Cn (cos(fi„r)|n, e) - i sin(0„r)|n + k, <?}) , (2.2) 

where r = gt is the scaled time and fi n is the generalized Rabi frequency |(| (in units of 
w ) 



0„ = y/[n + l)[n + 2)---{n + k). (2.3) 

In the large photon number case, namely when the average photon number n is large 
enough (fi> fc), the Poisson distribution {C n } is mainly concentrated near n and we have 
n ~ ft. 

We now specialize to the case k = 4. To facilitate an analytical treatment, we write 
the generalized Rabi frequency O n as 

2 / 10 35 50 24 , x 

n n = nWl - x, x= — + — + — + — . 2.4) 

n n z n 6 n 4 

For the large photon number case, a; is a small term and we can use the following Taylor 
expansion 

VT^=l + |-^ + .... (2.5) 
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Here we can only neglect the terms in f2„ which are much smaller than 1, namely, we 
must keep the terms to n~ 2 in (2.5). This is because fl n is the argument of trigonometric 



functions in the wave function (2.2) and a small constant (for example tt/2) can effect the 
physical quantity drastically. We therefore need to calculate the first three terms in ( |2.5| ) . 
We then obtain the approximate generalized Rabi frequency 

Q„ n 2 + 5n + 5, fl n -4 » n 2 - 3n + 1, (2.6) 

and fi„ — H„_4 = 4(2n + 1). Here the constant terms (5 or 1) are much smaller than n 2 
but cannot be neglected. 

3 Photon number distribution 

Taking the trace in the atom space, one finds for the reduced density operator of the cavity 
field 

OO 

Pf = ^ c nC* n [cos(tt n t) cos(Q m t)\n) (m\ + sin(fi„t) sin(fi m i)|n + 4) (to + 4|] , (3.1) 

m,n— 

from which the photon number distribution (PND) Pnij) = (w\pf\ti) is obtained as 

Pn{r) = |C„| 2 cos 2 (fi„T) + |C n _ 4 | 2 sin 2 (r! n _ 4 T) 

» \C n \ 2 cos 2 [fl n - 4 T + 4(2n + 1)t] + \C n -i\ 2 sin 2 [n„_ 4 r] . (3.2) 

It is easy to see that the PND is a 7r-periodic function of the scaled time r. This means 
that at r = 7r, we have P n {^) = |C„| 2 , a Poisson distribution. At t — tt/2, we have 
P„ (7r/2) = | C„ _ 4 1 2 , where we have used that fact that n 2 — 3n + 1 is always an odd 
integer. This means that the PND is also a Poisson distribution, but displaced by 4. 
However, the wave function (pure state) is not necessarily a coherent state; we shall see 
in the next section that the wave function is actually a Kerr state [[7j. 
When t = 7r/4, we have 4(2n + l)7r/4 = (2n + 1)tt and therefore 

P n (t) = |C„| 2 cos 2 [{n 2 -3n+ 1)tt/4] + |C„-4| 2 sin 2 [{n 2 - 3n + 1)tt/4] . (3.3) 

Writing n = 4m + i, i = 0, 1, 2, 3, we have 

(n 2 - 3n + 1) J = (4to 2 + 2im - 3to)tt + j £ s ^ J jj J' ^ (3.4) 



Then the PND is finally obtained as 



Pn{t) = \{\C n \ 2 + \C n ^\ 2 ), (3.5) 



namely, the average of PND at r = and 7r/2, as shown in Fig.l. 



When r = 7r/8, fi„ — f2„_4 = mr + ir/2 and the PND is obtained as 



P n (t) = (|C„| 2 + |C„- 4 | 2 ) sin 2 (n 2 - 3n + 1) 



7T 



(3.6) 
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which is a strongly oscillating function; in other words, the photon number distribution 
exhibits strong oscillation. Writing n = 8p + i, i = 0, 1, • • • , 7, we can further write ( |3.6| ) 



So for any n the photon number distribution is non-vanishing and, in that sense, the 
oscillation is not perfect. However, as can be seen from Fig.l, at the slightly earlier time 
r = 7r/8 — £, the photon number distribution exhibits perfect oscillation in that it becomes 
zero at n = 8p and 8p + 3. 

4 Entropy and pure states of the cavity field 

The entropy S of a quantum-mechanical system is a measure of how close the system is 
to a pure state and is defined by pi p) 



where p is the density operator of the quantum system and the Boltzmann constant is 
assumed to be unity. 5 = for a pure state and S > for a mixed state. In our model, 
the initial state is prepared in a pure state, so the whole atom-held system remains in 
a pure state at any time t > and its entropy is always zero. However, due to the 
entanglement of the atom and the cavity held at t > 0, both the atom and the held are 
generally in mixed states, although at certain times the held and the atomic subsystems 
are 'almost' in pure states. 

Since the initial state is a pure state, the entropy Sf of the cavity held equals the 
atomic entropy The entropy Sf or Sa, which is referred to as the entanglement of 

the total system in quantum information, is used to measure the amount of entanglement 
between the two subsystems. When Sf = Sa = 0, the system is disentangled or separable 
and both the held and atomic subsystems are in pure states. 

It is more convenient to calculate the entropy of the 2-level atomic system. From 



as 




(3.7) 



S = -Tr[pln(p)}, 
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Eq.(2.2) the atomic reduced density operator pa can be readily obtained as 

p A (T) = Tr F (\^(T))(^(T)\)=p 11 \g)(g\+p 12 \g)(e\+p 21 \e)(g\+p 22 \e)(e\, (4.8) 

where 

oo oo 

= \Cn\ 2 sm 2 (n n t), p 22 (r) = Y, \C n \ 2 cos 2 (n n t), 

oo 

Pi2(t) = P21 = Y C n+^ c n cos(il n+4 t) sin(n„t). (4.9) 
The field and atomic entropy Sa = — Tta(Po. hi(/0 a )) can be expressed as 

S F = S A = -tt+ ln(7r + ) - tt_ m(7T_), (4.10) 
where ir± are the eigenvalues of the atomic reduced field density operator pa 



7T± = -(I±y/{P22-Pll) 2 +Mpi2\ 2 )- (4.11) 

Noting that f2„ is an odd integer for large photon number, it is easy to see that the 
entropy Sf of the cavity field is a 7r-periodic function of r. We plot the entropy Sf as 
a function of r for n = 50 and < r < it in Fig. 2. We observe that the entropy is 
dynamically reduced to zero at it/ 2 and the cavity field can be periodically found in pure 
states. We also observe that the entropy falls quickly to minima near r — 7r/4, 3tt/4 (see 
also Fig. 3). 

Those calculated features can be explained analytically. The periodicity of Sf implies 
that Si?(m7r) = for any positive integer m. In this case the wave function ( |2.2D can be 
readily obtained as 

y(mir) w e lm7I \a) ® |e); (4.12) 

namely, the cavity field is in the coherent state \a) (up to a phase). 

When r = 7r/2, we have p\\ — 1, P22 = P12 = and thus the entropy Sf(jt/2) = 0. In 
this case the wave function of the system is 

*(tt/2) = -i (jr C n (-l) n(n+1)/2 \n)\ ® |.g) = e —/ 2 | - o.tt) ® |.g); (4.13) 

Vn=0 / 

that is, the atom is in its ground state and the field is in a pure state, which is known as 
a Kerr state 

00 n 

|a, 7 > = e-l Q l 2 / 2 Y " e *m(»-i)/2|„). (4.14) 



When the atom is detected in its ground state, the cavity field will be in the Kerr state. 

When r = 7r/4 and 37r/4, we have pn = p-ii and pyi ~ 0, where we have used the 
following fact for large photon number 



cosh(n) = 2_, 



EE. n 2n+1 



E 



sinh(n) 



(4.15) 



J2n)\ £j(2n+l)! 
The entropy is then Sf = — hi(2) ss 0.67315 (The exact numerical result is 0.69314, see 
Fig. 3). This means that the cavity field is in a mixed state. However, at a slightly earlier 
or later time r = tt/A + S r , the entropy is dramatically reduced to a minimum, as can be 
seen in Fig. 2 and Fig. 3. For the specific choice 



S r = — =rSi, r = ±l,±3,-- 
Ion 



(4.16) 



we obtain 



(O n - ft„_ 4 ) (J + <5r) M (2n + r + 1)tt + tt/2, (4.17) 
where we have used the relation n ~ n for large photon number h. In this case we can see 



that the entropy is approximately zero. In fact, at r = 7r/4 + 8 r , using (4.17) we obtain 
the following factorized state (writing a = \a\e l( ^) 



|*(tt/4 + *,.)) 



^C n sin(O n r)|n) ® (e 4 ^|e) + i\g)) 



= «5(f+5 r ) 



\n=0 



-i6<5 r 



)7 r/2 + 25 r ) - e - 25 (f+^) \iae- l6Sr ,-ir/2-2S r ) 



(e^\e)+i\g)) 



(4.18) 



from which we see that the cavity field is a macroscopic superposition of two Kerr states 



1 



V2 



,i5{f+S r ) \-i ae -iMr^/ 2 + 2S r ) - e~ i5 ^ +<5 ") \iae- m % -tt/2 - 2<5 r )) . (4.19) 



The photon number distribution can be derived from the above wave function. Here 



we write down the PND directly from ( p.l[ ) and ( 4.17 ) 

P n (ir/4 + d r ) = (|C„| 2 + |C„-4| 2 )sin 2 [{n 2 - 3n + 1)(tt/4 + S r )] . (4.20) 

In Fig. 4 we compare this approximate result with the numerical result. We see that the 
agreement is good. 



We now show that the entropy is approximately vanishing. Using \C n 



|CL_ 



4 m 



P^0| ) we find 

oo oo 

1 = p n(K/4: + S r ) = 2^2 \C n \ 2 sin 2 [(n 2 + 5n + 5)(tt/4 + S r )] = 2p u , (4.21) 



n=Q 



from which we have pn(7r/4 + S r ) = 1 — /?22(7r/4 + ^r) = 1/2. ^From (4.17) we have 

pu (n/ 4 + 6 r ) = -e 41 Vi (tt/4 + <5 r ) = -e 4l */2. (4.22) 
So the entropy 5^(^/4 + <5 r ) w 0, as expected. 

5 The Q-function and Population Inversion 

The quasi-probability distribution Q-function is defined as |l(| : 



Q{0) = -(P\Pf 



(5.1) 



where pf is the reduced density operator of the cavity field given in Eq. (3.1) and \j3) is 



a coherent state. Inserting Eq. (3.1) into Eq. (5.1) we can easily obtain the Q-function of 
the cavity field 



QW) 



\p\ 



—j=Cn cos(fi„r) 
n =o v n! 



E // , jN| CnCOs(» tt r) 



(5.2) 



When r = and tt, the cavity field is in a coherent state and its Q-function is Pois- 
sonian, as can be seen from Fig. 5. At r = 7r/2, the Q-function has two components. At 
t = 7r/8, the Q-function has 8 components, and the field is in a mixed state (its entropy 
is 0.6888). So the interference between components results in imperfect oscillation of the 
PND. 

When r — n/4: the Q-function is composed of 4 well-separated components, the entropy 
reaches its maximum and the field is in a mixed state. In a short time r = 7r/4+5i, the field 
is nearly pure and the Q-function separates into 8 components, the interference between 
which resulting in strong oscillation of the PND, as can be seen from Fig. (4). However, 
at r = 7r/4 + ^3 and r = tt/4 + 65 the cavity field is less pure and the Q-function has less 
well-separated components, and the oscillation of the PND is less marked (see Fig. 4). 

Finally, we briefly describe the atomic population inversion (API) , given by 



W(t) = (o- 3 ) - J2 \Cn\ 2 cos{2il n T). 



(5.3) 



n=0 



In the one or two photon cases this exhibits the usual periodic collapse and revival phe- 
nomenon. The 4-photon system also shows this behaviour. In Fig. 6 we give a plot of the 
API (03) for n = 50, from which one can see that the API for the 4-photon system also 
exhibits periodic collapse and revival. However, due to the nonlinearity of the system, this 
effect does not last as long as in the one or two photon cases. 



6 Conclusions 



In this paper we studied the four-photon JC model and its dynamical behaviour. In the 
large photon number approximation, the Rabi frequency fl n depends quadratically on n 
and this leads to nonlinear effects in the cavity. 

By examining the entropy of the field or atom we found that the system is approxi- 
mately disentangled at certain times. In this case both the atom and cavity field are in 
a pure state and those pure states are explicitly given. It is interesting that the cavity 
field exists in both a nonlinear Kerr state and a macroscopic superposition of Kerr states, 
analogous to the coherent state and the Schrodinger cat states in the one or two photon 
JC model. One may refer to this analogous Kerr superposition as a Kerr cat state. 

We analysed the PND in the large photon number approximation. At certain times 
the PND exhibits strong oscillation. From the Q-function we know that the cavity field 
at those times has eight components; it is the interference between those components 
which causes the strong oscillation of the PND. Note that these oscillations depend more 
sensitively on the interaction time than in the one or two photon cases due to the quadratic 
dependence of the generalized Rabi frequency on photon number in the present work. 

We also considered atomic number inversion and found that the phenomenon of peri- 
odic collapse and revival occurs; it is however a short-lived phenomenon due to the effects 
of nonlinearity. 

Experimental implementation of the theoretical idea herein proposed would possibly 
depend on the use of a trapped neutral atom or ion in a high-Q cavity, or an atomic 
beam in transit through the cavity |DJ. In any realistic scheme off- resonant effects must 
be taken into account, and for this reason the atomic species used must provide a good 
approximation to a two-level system. 
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Figure 1: Photon number distributions for n = 50 at r = 0,7r/8,7r/4 and 7r/24. For the 
r = 7r/8 case, the photon number distribution at slightly earlier time r = ir/8 — 7r/24000 
(dashed line) is also presented. 
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Figure 2: Entropy Sf{t) of the cavity field against the scaled time r for n = 50. 
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Figure 3: Entropy Sp(r) of the cavity field around t — tt/A for n = 50. The vertical lines 
with labels r — ±1, ±3, ±5 represent r = 7r/4 + rSi, Si — 7r/800. 





Figure 4: Photon number distribution of the cavity field at r = ir/4 + Si and r = 7r/4 + S$ 
for n = 50. The solid line is the exact numerical result and the dashed line is from the 
approximate result Eq.(4.2C). 
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